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■ Abstract 

The statics and dynamics of a surface separating two phases of a relativistic 
quantum field theory at or near the critical temperature typically make use of a 
free energy as a functional of an order parameter. This free energy functional also 
affords an economical description of states away from equilibrium. The similarities 
and differences between using a scalar field as the order parameter versus the energy 
density are examined, and a peculiarity is noted. We also point out several concep- 
tual errors in the literature dealing with the dynamical prefactor in the nucleation 
rate. 
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1 Introduction 



Different phases of matter are separated in space and time by dividing layers called do- 
main walls. The dynamics of the domain walls is the determining mechanism governing 
phase transitions in an intermediate range between very slow, quasi-static and very rapid, 
dynamical processes. In this intermediate range the phase transition speed and the speed 
of external constraints are comparable to each other. 

The dynamics of phase transitions is an involved subject even in macroscopic sys- 
tems. First of all, phase transitions can be different. They may include slow burning 
or deflagration, detonation, condensation, evaporation, and many other forms of tran- 
sition. Nevertheless, the basic conditions of all these transitions have some similarities. 
These arise from the basic conservation laws and from the requirement of local, or at least 
approximately local, equilibrium. 

In a dynamical situation the approach using the equation of state with a first order 
phase transition is identical both in compression and in expansion. If the compression is 
supersonic, shock or detonation waves are formed where the final new phase is immediately 
created. The phase transition speed influences only the width of the shock front. On the 
other hand, for slow dynamics and rapid phase transition the shock front width is primarily 
determined by the transport coefficients, viscosity and heat conductivity, and not by the 
phase transition speed. 

At high energies a relativistic treatment is frequently necessary. It is important to 
mention that a system is also relativistic if the matter is radiation dominated, meaning 
that the rest mass of the constituent particles is zero or negligible compared to the fourth 
root of the energy density. These types of systems must be treated as relativistic even if 
the collective velocities are small. This was one of the important new features recognized 
in ref. [0. If conserved charges do not exist the application of the conventional theory of 
phase transition dynamics []2], |3|] is not possible, and this leads to essential differences in 
the phase transition dynamics. Most importantly the flow is tied to the energy flow; thus 
energy or heat current and the fluid flow are identical and heat conduction (with respect 
to the flow) may not occur. Therefore the coefficient of heat conductivity does not exist, 
and cannot govern the phase transition speed. 

The spatial configuration of phase transition instabilities varies in large homogeneous 
systems. When the conditions for the occurrence of a new phase are established the most 
common form of the appearance of the new phase happens via the formation of small 
critical size bubbles or droplets. (Sub-critical size bubbles or droplets will shrink and 
vanish but supra-critical size ones will grow.) This configuration is called homogeneous 
nucleation. 

As the level of supercooling increases and the amount of the new phase increases 
other geometries become energetically more favorable, such as elongated cylindrical ob- 
jects (spagetti) or layers (lasagna). The phase transition dynamics is then referred to 
as spinodal decomposition, indicating that systems which supercool (or superheat) and 
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reach the adiabatic or isothermal spinodals on the phase diagram start the formation of 
the new phase in these configurations immediately. 

Finally, when the two phases are about equally abundant and/or the transition is 
extremely rapid the two phases form a somewhat random occupation of the configuration 
space called percolation. 

Most frequently explicit dynamical calculations are performed for the homogeneous 
nucleation geometry as this is usually the initial and the slowest of all. Nevertheless, the 
domain walls and their dynamical properties play an important role in all of the above 
mentioned configurations. One important feature of the nucleation studies that is not 
always recognized is that within the dynamical domain wall or droplet wall we do not 
have thermal equilibrium. Here, quantities like energy, entropy or particle density make 
sense, while other quantities that rely on thermal equilibrium, such as the equation of 
state, may not. 

If the fluid is perfect and consists of a single phase of matter in full phase equilibrium, 
the flow is adiabatic ||]. Nevertheless entropy production is still possible if the system is 
out of phase equilibrium @], or out of thermal equilibrium. Thus, in domain walls entropy 
may be produced even in perfect fluids as here we do not have an equation of state which 
would satisfy the requirements of equilibrium thermodynamics. 

In this paper we would like to discuss connections between the two most frequently 
encountered approaches to phase transitions using effective field theoretical models and 
phenomenological thermal and fluid dynamical methods. In section 2 we first demonstrate 
with a simple effective field theoretical model how one can obtain a Landau-model type 
of phenomenological description of critical dynamics. In section 3 we elaborate on the 
assumptions and error of different approaches used recently in the literature. We present 
a summary and conclusion in section 4. 

2 The Free Energy Functional 

A common approach to dealing with a phase transition in relativistic quantum field theory 
is to assume a uniform condensate (ft for some scalar field and then to compute quantum 
and thermal fluctuations about that condensate ||. This is especially so for systems 
exhibiting spontaneous symmetry breaking where the condensate field serves as an order 
parameter to distinquish the two phases. Examples include the Higgs field in electroweak 
theory and the sigma field or scalar quark condensate in strong interaction physics. The 
resulting free energy density / is obtained from the partition function Z in the usual way. 

f(<f>,T) = -T\nZ((f>,T,V)/V (1) 

Here V is the volume. Since the high temperature symmetric phase usually corresponds 
to = 0, it is customary to define the effective potential as the deviation of the free 
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energy from its value in the symmetric phase. 

V&T) = f(<f>, T)-/(0,T) (2) 

The point = is either a global minimum or only a local minimum depending on 
whether T is greater than or less than the critical temperature T c , respectively. Therefore 
the results of computations are often found to be, or parametrized as, a fourth order 
polynomial in 0. 

V(<j>,T) = J2 a n(TW (3) 

n=2 

Occasionally one will find additional terms of order 4 In (4>/A) or higher powers of 0, but 
these are not common and will not change our analysis qualitatively. 
To give a specific example we will use the parametrization of ref . . 

Vfa T) = i 7 (T 2 - T o 2 )0 2 - i«T0 3 + (4) 

The 7, T 2 , a and A are temperature independent constants, to be specified shortly. In an 
equilibrium state the free energy density is the negative of the pressure, f(T) = —p(T). 
Normalizing to the high temperature symmetric phase, and denoting the equilibrium 
pressure of that phase by Ph(T), we have 

/(0,T) = -p fc (T) + V(0,T). (5) 

It is not enough to specify the effective potential; the equilibrium pressure as a function 
of temperature must be specified too. 

If one knows the free energy as a function of T then the energy density e can be 
calculated from the basic thermodynamic identities. In particular, if we want to know 
the energy density for a specific value of the condensate field, including out of equilibrium 
configurations also, then we must hold it fixed during the temperature differentiation. 

e(0,T) = /(^T)-T^,T) 

= e,(T)-i 7 (T 2 + T o 2 )0 2 + ^ 4 (6) 

The first line above is the general thermodynamic identity; the second line applies to the 
specific parametrization under discussion. If one is at a local minimum of the effective 
potential, either (fih(T) = or <f>i(T) > corresponding to the high and low temperature 
phases, respectively, then it doesn't matter whether the order parameter is held fixed 
during the differentiation or not because df(<j),T)/d<f) = at those points. These two 
points correspond to thermodynamical equilibrium. 

If the system undergoes a first order phase transition at a critical temperature T c then 
at that temperature / will have two degenerate minima. For a finite range of temperature 
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above T c there will persist a higher, local minimum at (f>i(T) > 0. For a finite range of 
temperature below T c there will persist a higher, local minimum at <ph{T) = 0. The 
most important physical quantities are the latent heat L = Ae(T c ) = e h (T c ) — ej(T c ), the 
correlation lengths £~ 2 = d 2 f (0 , T c ) / d<f) 2 = d 2 f((j)i(T c ),T c )/d<f) 2 , which are equal for the 
fourth order effective potential at T c , and the surface energy a. At the critical temperature 
the planar interfacial profile has a nice analytical solution on account of the fact that the 
effective potential becomes symmetric. In the usual way one finds that the profile field 
<f>(x) is 

(x) = ^1 [i _ tanh(x/20] • (7) 

This interpolates through intermediate non-equilibrium states from </>/(T c ) for x <C — £ to 
4>h(T c ) = for x ^> £. The surface energy can then be expressed as 

,-r\ 2 



Xoo i^dx) 




The parameters in the effective potential can be expressed in terms of these physical 
parameters as follows. 

7 = 

rp2 

1 O 

a = 

For purposes of illustration, and with the QCD phase transition in mind, we choose the 
following numerical values. 

^ = T c 

L - ^ (10) 

For the pressure in the high temperature phase we take 37 effective massless bosonic 
degrees of freedom corresponding to gluons and two flavors of quark, plus a bag constant 
B — L/A to simulate confinement. 

Tin 2 

Ph (T) = ^-T*-B (11) 
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The numerical value for the latent heat L corresponds to a transition from the 37 degrees 
of freedom mentioned above to 3 degrees of freedom for massless pions. It is convenient to 
make everything dimensionless by measuring energies in units of T c and lengths in units 
of 1/T C , which practice we adhere to in the rest of this section. 

The effective potential (4) as a function of is plotted in figure 1 for three different 
temperatures: the critical one and 1% above and below it. This displays the traditional 
behavior of a strong first order phase transition. The plot includes negative values of 
which may or may not be allowable, depending on the origin of and its physical 
interpretation. 

The free energy density (5) as a function of is plotted in figure 2 for the same 
temperatures. The only difference between this and figure 1 is the addition of the pressure 
as evaluated in the high temperature symmetric phase. 

The energy density (6) as a function of is plotted in figure 3. It has a maximum at 
= 0. This follows directly from the expression for e given earlier: For small values of 
the deviation from e^(T) is quadratic in with negative curvature. Note that the energy 
density goes negative when is greater than about 2.0. 

Rather than using the field as the order parameter one might consider other choices 
more appropriate for the problem at hand. For example, one could take as the order 
parameter the deviation of the energy density from its equilibrium value and expand 
the free energy in a power series in this difference. This is the Landau approach to 
the description of fluctuations and to departures from equilibrium states Such an 
approach to the nuclear liquid-gas phase transition was implemented by Goodman et 
al. |§. It is also this same function that was used in working out the nucleation of 
relativistic first order phase transitions by two of us ffl. With these motivations we plot 
the free energy versus the energy in figure 4 by using as a parameter. There are a 
number of points to be made concerning this figure. First, there are two minima at 
each of the chosen temperatures. The minima occur at the energy densities e/(T) and 
e/ l (T). The lower of these two represents the equilibrium state while the other represents 
a metastable state except at T c when they are degenerate. The value of / at a minimum 
is equal to the negative of the pressure in that particular phase: Pi(T) = —f(ei(T),T) and 
Ph{T) = —f(eh(T),T). Second, the value of e;(0.99T c ) < 0! This is an indication of the 
inadequacy of the specific parametrization of the effective potential at this temperature, 
and probably at lower temperatures too. This is an important point to always be aware 
of when writing down any formula for an effective potential. Third, the free energy has 
a cusp at the location of the high density symmetric phase. This is an unavoidable 
consequence of the inference of the free energy from the effective potential. From figure 
3 or from the corresponding formula it is clear that the energy density is a maximum at 
= 0. This means that the plot of the free energy ends at the equilibrium density &h{T) 
of the symmetric phase with a cusp. Physically this must be a restriction of fluctuations 
to those associated with and only 0. This is too restrictive by far; there are certainly 
other physical processes not taken into account. For example, placing particles in a box 
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in contact with a heat reservoir allows for the exchange of energy between the particles 
and the reservoir, resulting in fluctuations in energy and pressure. These processes are 
always present in systems at fixed temperature and are not accounted for, or associated 
with, fluctuations in an order parameter. 

Now let us follow the Landau construction of the free energy away from equilibrium 
states using the energy density as the order parameter with no reference to </> whatsoever. 
This is the construction made in |i| in the context of constructing spherical surfaces 
separating two phases at noncritical temperatures. The total free energy of interaction is 
the sum of a gradient energy and a free energy density integrated over space. 



Fj{e(x)} = J d 3 x 



^(Ve(x)) 2 + /(e(x),T) 



(12) 



The / is most economically written as a fourth order polynomial in e. 



f(e,T) = J2 b n(T)e n (13) 

n=0 

The coefficients are functions of temperature. They are determined by several require- 
ments. The first is that / has minima located at ei(T) and eh{T), such that / at those 
points be equal to the negative of the corresponding equilibrium pressure. 

f(e h (T),T) = - Ph (T) 

f{e l {T),T) = - Pl (T) (14) 

This results in four equations involving the five b n . Two more equations result from 
fixing the correlation length (which is the same in the two phases only at T c ) and the 
surface energy o (which, strictly speaking, is only well defined at T c ). This set of six 
equations is not inconsistent because the coefficient of the gradient energy, K, must also 
be determined. 



d'f/de% =eh(Tc) = d'f/de% =ei(Tc) a2 



K = JY ( 15 ) 

We insist upon the same equilibrium energy densities and pressures as used in the effective 
potential approach, and use the same correlation length and surface energy at T c too. 

e(x) = -[e h (T c )+e,(r c ) + Ae(T c )taoh(a;/20] (16) 




a = K dx\ — \ (17) 
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Then the Landau expansion of the free energy as a function of the energy density is 
obtained and is plotted in figure 5. The f(e, T) is now a smooth and well-behaved function. 
Its only failing is that the equilibrium energy density at 0.99T C is negative, but that is 
a consequence of insisting that it be the same as for the effective potential, which was 
negative. (Actually, fluctuations into states with negative energy can be avoided if one 
uses the Laurent expansion instead of the fourth order polynomial approximation |J, 
but only if the lower minimum is at positive energy density.) In fact, with the Landau 
approach we have a closer connection between the physical observables and the free energy. 
For example, we can easily specify the equilibrium energy density and pressure in each 
phase, making it simple to avoid such unwelcome behavior as a negative energy density 
in equilibrium. Of course the direct connection with the condensate field is lost. 

The main difference between the two approaches is that the field theoretical approach 
parametrizes the out-of-equilibrium configurations in terms of an effective potential as a 
function of <fi that shows smooth, quadratic minima around the two equilibrium states, 
whereas the Landau type of thermodynamical approach does the same in terms of the free 
energy as a function of the energy density e. These two approaches are not fully equivalent: 
the highly nonlinear /(e) dependence obtained from this particular field theoretical model 
does not yield a smooth quadratic minimum in terms of e at (where = 0). This leads 
to unphysical estimates for energy density fluctuations, in our opinion. Which approach 
one takes depends on the physical situation being addressed. 

3 Dynamics of First Order Phase Transitions 

A phase transition occurs because of some change in the global properties of the system. 
Examples of such changes are: an expansion chamber or a sudden quench in temperature 
in the laboratory, the expansion of hot matter in a high-energy heavy ion collision, and the 
expansion of the early universe. Generally we need to compare the phase transition rate to 
the rate of expansion or quench. The speed at which the phase transition proceeds comes 
into play when the speed at which the external variables, such as the volume, becomes 
comparable to or exceeds that of phase conversion. 

If the rate of change of an external variable is slow there is sufficient time to maintain 
phase equilibrium on a coarse-grained time scale. This also means that all other equili- 
bration processes are completed as these invariably require less time and less interaction 
than the conversion of one phase of matter to another. 

Thus, in the case of slow external dynamics and rapid phase equilibration the matter 
is in complete equilibrium, including phase equilibrium, and the equation of state of the 
matter undergoing a first order phase transition is given by the Maxwell construction. 
Then we have a fully developed mixed phase, and the phase abundances are typically 
given by the conservation of energy and entropy. No information on dynamical processes 
is needed. This type of transition is adiabatic. Even in moderately fast dynamical ex- 
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pansions there are only small deviations from the ideal and complete phase equilibrium 
(Maxwell construction). This deviation results in some delay in the creation of the new 
phase, leading to supercooling or superheating and extra entropy production. [Q 

For heavy ion reactions the first attempt to explicitly evaluate the phase transition 
speed of homogeneous nucleation was described in fH, pTOfl . The homogeneous nucleation 
mechanism correctly describes the initial stage of the phase transition when the abundance 
of the newly created phase is still small and when the rate of phase conversion is the slowest 
process. 

Here a couple of remarks are necessary. There are several requirements to form critical- 
sized bubbles or droplets of the new phase. Both pressure balance and temperature 
equality should be established between the phases; this requires the transfer of energy 
and momentum across the phase boundary. If local equilibrium is assumed both before 
and after the formation of the new phase then we cannot relax the requirement of pressure 
and temperature equilibrium. 

Langer's modern theory of nucleation yields the following formula for the rate: 

/ = f fioe~ AF/T (18) 

where AF is the change in the free energy of the system due to the formation of the critical 
droplet. Qq is a statistical prefactor which measures the available phase space volume, k 
is a dynamical prefactor which determines the exponential growth rate of critical droplets 
which are perturbed from their quasi-equilibrium radius R*. 

k = j t HR(t) -R*} (19) 

The dynamical prefactor has been calculated by Langer and Turski [0, and by Kawasaki 
TT|] for a liquid-gas phase transition near the critical point, where the gas is not dilute, 



to be 

K =ftPT (20) 

This involves the thermal conductivity A, the surface free energy a, the latent heat per 
molecule i and the density of molecules in the liquid phase ri£. The interesting physics 
in this expression is the appearance of the thermal conductivity. In order for the droplet 
to grow beyond the critical size, latent heat must be conducted away from the surface 
into the gas. For a relativistic system of particles or quantum fields which has no net 
conserved charge, such as baryon number, the thermal conductivity vanishes. The reason 
is that there is no rest frame defined by the baryon density to refer to heat transport. 
Hence this formula obviously cannot be applied to such systems. The dynamical prefactor 
for such systems was first evaluated a' la Langer by two of us to be 



4(7 

K 



(Aw) 2 Rf 



4 

3* + C 



(21) 
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Here r\ and ( are the shear and bulk viscosities, respectively, and Aw is the enthalpy 
density difference between the two phases (equal to the latent heat since the pressures 
are equal at the critical temperature). This fully relativistic expression was subsequently 
generalized by Venugopalan and Vischer El to include a conserved baryon number. 
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K 



XT + 2 ( -7J + C 



(22) 



In the non-relativistic limit and when the viscosities are small compared to the heat 
conductivity this reduces to Langer's expression. In the relativistic limit with no net 
baryon number, where effectively A — ► 0, this reduces to equation (21). 

Unfortunately there are several erroneous expressions for k which subsequently ap- 
peared in the literature. In [13 one finds that in the relativistic limit with no net baryon 
number 



K = ^j^m ■ (23) 

This was "derived" under the assumption of vanishing shear viscosities. An expression 
for k in the nonrelativistic limit with non-zero baryon number was found which differed 



from Langer's result too. By including the viscosities another paper JTj] finds that 



2awh 



(Aw)mi 



+ 



a 



c 2 s (Aw) 2 Rl 



(24) 



Here c s is the speed of sound in the low density phase 
vanish this reduces to the result of [l - 
to the result of PJ albeit with a factor of 1/c 2 rather than 4 



In the limit that the viscosities 
and in the limit that they are large this reduces 



The results found by |T3[ and [14] are suspicious because they predict a dynamical 
growth factor even in the absence of viscous forces. These results are wrong because of a 
misinterpretation of a derivative. Both papers follow the approach of |l[ relatively closely 
except for one crucial point. At some point in the analysis one encounters the derivative 
df(e,T)/de evaluated at the radial profile e = e(r) of the energy density for a spherical 
bubble or droplet connecting the two phase. Thus e(r) varies between q and e^. The 
free energy density varies accordingly. As explained carefully and in some detail in PJ, 
this derivative is taken at fixed temperature. At the equilibrium points it is zero: see the 
discussion in section 2 above. However, both |H| and JHJ equate / with — p irrespective 
of whether we are at one of the minima or not. They then assume an equation of state 
p = c 2 e to finally obtain df(e,T)/de = —c 2 s . From that point on the analysis diverges, 
and the erroneous results for k inevitably follow. 

It is also possible to argue that the erroneous results for k are not only mathematically 
incorrect but physically incorrect too. What is relevant is the slowest required process, 
not the fastest. Furthermore it is incorrect to add up the rates of all processes to obtain 
the highest possible phase transition speed. If not all required processes are completed the 
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phase transition is not complete either. For example to establish pressure (momentum) 
equilibrium O, [Ml is not sufficient, because then a subsequent step of establishing thermal 



equilibrium is required, and only when both processes are over has the phase transition 
been completed. 

In ref. [0J the slowest process, heat transfer via viscosity, was evaluated. In some cases, 
when we have non-negligible net baryon charge in our system, the heat conductivity may 
also contribute to the heat transfer to the new phase, thus speeding up somewhat this 
slowest of all processes [O. If the transport coefficients are all vanishing so that thermal 



balance cannot be achieved, we will never reach both phase and thermal equilibrium, so 
the rate tends to zero as stated in ref. [|TJ. 

In case other degrees of freedom exist which permit energy transport between the 

TBI, these processes can be combined with 



phases leading to a common temperature 
other transport processes like viscosity and/or heat conduction leading to faster tempera- 
ture balance and a higher rate. However, sound waves do not lead to dissipation as perfect 
fluid dynamics is adiabatic (in the absence of shock waves). Thus, the sound wave as the 
sole mechanism cannot characterize dissipative transport processes, as was mistakenly 
claimed in 
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A final remark: If the external dynamical evolution is much faster than the processes 
involved in the phase and kinetic equilibration then the matter involved in the phase 
transition loses both phase and thermal equilibrium, and we have to abandon the thermal 
and fluid dynamical approaches altogether. An example in the context of high energy 
heavy ion collisions is presented in [17] where an effective field theoretical approach is 



used. Generally this will be the case for small, rapidly developing systems only. 



4 Conclusions 

In this paper we have examined two subtle issues in the dynamics of domain walls. One of 
them is the difference between the effective potential as a function of a scalar field and the 
free energy as a function of energy density, both evaluated near the critical temperature. 
The other involves a misconception in the literature about the derivation of the dynamical 
prefactor, or growth rate, in homogeneous nucleation theory. It is important to understand 
these subtleties because of their importance in such physical environments as cosmology, 
astrophysics and high energy nuclear collisions. 
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Figure 1: The effective potential as a function of the field. 
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